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Abstract 

We consider a diffusion-limited reaction in case the reacting entities are not available simul- 
taneously. Due to the fact that the reaction takes place after a spatiotemporal accumulation of 
reactants, the underlying rate equation has to be modified by additional non-local terms. Owing to 
the delay effects a finite amount of reactants remains localized, preventing a further reaction and 
the asymptotic decay is terminated at a finite density. The resulting inhomogeneous non-zero sta- 
tionary concentration is stable against long wave length fluctuations. Below a critical wave vector 
k c the system becomes inert, whereas a complete decay is realized above k c . The phase diagram 
for the one species-annihilation process A + A — > exhibits a behavior comparable to a second 
order phase transition. Obviously the memory effects are equivalent to long range interaction and 
the non-local kinetics is basically independent on space dimensions. 
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I. INTRODUCTION 



In case the typical transport time td of reacting entities is much larger than the re- 
action time tr, the reaction process is affected by spatial correlations. Consequently the 
global concentration c(t) has to be replaced by a local density function c(r, t). When the 
reaction occurs within a cell, centered around the position r, the classical mean-field rate 
evolution equation is completed by a diffusive term, which disposes of low dimensional fluc- 
tuation corrections. The influence of such fluctuations on the longtime behavior had been 
demonstrated in a series of papers Q], for a recent review see Q. The newly emerging 
nonequilibrium behavior is extracted using more sophisticated methods as field theoretic 
treatment or extensive numerical simulations ^J. In the one species-annihilation process 
a particle A is annihilated upon encounter according to the reaction scheme A + A — > 0. 
The reaction, realized with a certain rate u, takes place immediately. The time evolution of 
the local concentration c(r, t) satisfies the rate equation 

d -^ = DW 2 c{v,t)-uc\r,t) (I) 

Suppose that I is a typical length scale until reactants meet and c ~ l~ d a typical concen- 
tration in that region, the condition td ^> tr leads to 

l 2 - d » - (2) 
u 

Since the diffusion constant D and the reaction rate u become fixed by microscopic processes, 
the spatial variation of the concentration on the length scale I should be relevant for d < 
2. Due to the short reaction time tr the particles annihilate very rapidly. Therefore a 
subsequent reaction can only take place in turn the particles has passed through the distance 
/. This effect slows the density decay down. Asymptotically one finds c(t) ~ t~ d l 2 for d < 2 
and logarithmic corrections at the critical dimension rf c = 2 p. Recently, the result has 
been generalized to multispecies pair annihilation leading to a modified exponent for the 
asymptotic decay In the opposite case d > 2, spatial fluctuations are irrelevant and the 
mean-field equation provides a valid description. The diffusion-limited reaction Eq. ([TJ is 
based on the assumption that the reactants are available simultaneously. In the present paper 
we emphasize, that the generic behavior of the system may be changed when additional non- 
linear delay effects are included into the consideration. In particular, we discuss a reactions 
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which takes place after a sufficient accumulation time and after cumulating reactants within 
a spatial region. Under that conditions the reaction time scale tr is modified and likewise 
Eq. (J2J) is changed. Accordingly we demonstrate, the long time behavior is dominated 
apparently by the delay effects. The system is able to reach a stationary state with a finite 
concentration instead of exhibiting an algebraic decay in time as it follows using Eq. (JTJ. 
Although the crucial factors governing the dynamics of systems, comprising many "units", 
consists of interaction and competition, there is an increasing interest to include memory 
effects as a further unifying feature of complex physical [?], 0] as well as biological systems 
0. Our model can be grouped into the effort to discuss delay and feedback mechanism. It is 
well-known that evolution equations with memory kernels can be derived following the well 
established projector formalism due to see also jllj]. That approach had been applied 
successfully for the density-density correlation function studying the freezing processes in 
undercooled liquids Q] . Recently a Fokker-Planck equation with a non-linear memory 



term was used to discuss anomalous diffusion in disordered systems |14j|. The results could 
be confirmed by numerical simulations including diffusion on fractals see also lfil \\\. 



Moreover, it was argued |Tj| that mobile particles remain localized due to the feedback- 
coupling. Notice that a formal solution of whole class of non-Markovian Fokker-Planck 
equations can be expressed through the solution of the Markovian equation with the same 
Fokker-Planck operator Q|. The non-Gaussian fluctuations of the asset price can be also 
traced back to memory effects 20]. An additional cumulative feedback coupling within the 
Lotka-Volterra model, which may stem from mutations of the species or a climate changing, 
leads to a complete different behavior [2]] compared to the conventional model. If the 
Ginzburg-Landau model for the time evolution of an order parameter is supplemented by 
a competing memory term, the asymptotic behavior and the phase diagram is completely 



dominated by such a term 



22j. Whereas the feature of the approach, proposed in those 



papers, consists of self-organization, i. e. the time scale of the memory is determined by 
the relevant variable itself, for instance the concentration, there is a broad class of models 
with external delay effects 

That case is characterized by a given external memory kernel. The differences of both 



approaches will be discussed also in |27[ . The spreading of an agent in a medium with long- 
time memory, which can model epidemics, is studied in Time-delayed feedback control 
is an efficient method for stabilizing unstable periodic orbits of chaotic systems j^J. Time 
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delay may induce various patterns including travelling rolls, spirals and other patterns 



The influence of a global feedback is studied recently in a bistable system 31], where the 
purpose of that paper is a discussion of the domain-size control by a feedback. Even in an 
open quantum system non-Markovian dynamics is characterized by time-non-locality in the 
equation of motion for the reduced density operator [3^ . 

In view of the large variety of systems with feedback couplings it seems to be worth to study 
simple models, which still conserve the crucial dynamical features of evolution models as 
non-linearities and moreover, as a new ingredient, delayed feedback-coupling. In the present 
paper we discuss the influence of a non-Markovian memory term on chemical reactions. The 
retardation effects are characterized by memory kernels, which are chosen in such a manner, 
that they compete with the conventional diffusion and the instantaneous non-linear reaction 
terms. 



II. MODEL 

Although a chemical reaction is characterized by non-linear reaction terms, as it is in- 
dicated in Eq. the reaction process may be further modified when the reaction takes 
place after an accumulation time. In that case the time evolution of the concentration could 
also depend on the history of the sample to which it belongs, i. e. the changing rate of the 
concentration should be influenced by the changing rate in the past. Thus the evolution 
for global concentration c(t) has to be supplemented by memory terms. Such term model, 
for instance the way on which a seed concentration had been accumulated by a delayed 
transport mechanism, originated by the environment of the reactants. With other words, 
the changing rate of a certain quantity at time t is also determined by the accumulation 
rate at a former time t' < t. In between, i.e. within the interval r = t — t', the reactants 
are enriched while changing the concentration at t'. Regardless that process and further 
fluctuations the available amount of concentration at time t is governed by instantaneous 
loss term as well as on the changing rate at former times t' . Consequently the evolution 
Eq. (JTJ) should be modified according to 

t 00 

d t c{r, t) = V 2 c(r, t) - uc 2 (r, t) - J dt' J dr'K(r -r',t- t'; c) (3) 

-00 



4 



As the main ingredient we assume that the kernel K is determined by the concentration 
c and its derivative. With other words, the time and the spatial scale of the memory is 
governed by the spatiotemperal scale of the concentration. Under that conditions the chem- 
ical reaction becomes obviously a many-body problem. The reaction species are embedded 
into an environment of all the other particles of the system which give rise to feedback and 
memory effects. 

Let us now specify the memory kernel introduced in Eq. Q by the expression 



l oo 

d t c(r, t) = DV 2 c(r, t) - fi J dt' J dr'V 



c(r-r'A-t') ■dt.VdT't' 



-oo 
t oo 



- uc 2 (r, t) — X J dt' J dr'c(r -r',t- t')d t ,c{r', t') (4) 



o 



The second term on the right hand side characterizes a delayed diffusion with diffusion 
parameter [i. This term conflicts to the conventional local diffusive one with the parameter 
D. Large concentration fluctuations at both times t and t' contribute to the behavior of the 
system. Different to the constant D the sign of the parameter [i can range between positve 
and negative values. Likewise the third and the fourth term exhibit an extra competition 
between the local and the accumulated as well as the delayed reaction, where we assume 
u > and A arbitrary posisitive and negative. Using scaling analysis it is easy to verify that 
the last term with the coupling parameter A is relevant in all dimensions. In contrast the 
spatial dependence of the local w-term is only relevant for d < 2. The diffusion term and 
the delayed diffusion term are of the same order of magnitude. Follwing the line resulting in 
Eq. (j2J), we can estimate the characteristic time scales. Whereas the diffusion time remains 
unchanged, To ~ l 2 /D, there occurs an effective reaction time tr which is composed of the 
parameters /i, u and A of the three non-linear processes, introduced in Eq. (@J). The effective 
reaction time is fixed by 

1 



l d \ + l d ~ 2 fi 
The condition td 3> tr yields 



u 



+ + l d (l d \ + fil d ~ 2 ) (5) 



A/ 4 + fil 2 + l 2 ~ d Du > D 2 . (6) 

In the special case A = \x = this relation leads anew to Eq. (J2J). For non-zero memory 
parameters [i and A Eq. © is fulfilled in each space dimension. Insofar, the memory mimics 



a kind of very long-range forces. In the forthcoming section, we search for the stationary 
solutions in arbitrary dimensions d. 



III. STATIONARY CONCENTRATION 

In this section we show that the asymptotic behavior of the concentration is changed 
significantly when the retardation effects are taken into account. Firstly let us consider the 
case of vanishing memory parameters ji = and A = 0. In mean-field approximation the 
density decay obeys c ~ 1/ut. For d < 2 the density fluctuations leads to an additional 
time dependence of the parameter u — > u(t). The explicit time dependence of u can be 
estimated easily by scaling arguments. The parameter u scales as u ~ l d ~ z where the 
dynamical exponent z is related to the time scale by I ~ t x l z . Because for conventional 
diffusion the dynamical exponent is z = 2 we find u(t) ~ t d /2-i which leads immediately to 
the asymptotic relation c(t) ~ t~ d l <1 in accordance with more refined methods 0, compare 
also |2J . All particles are annihilated due to the reaction. This picture is changed, when the 
delay effects are enclosed into the consideration. The asymptotic behavior in time is obtained 
by applying Laplace transformation. Moreover let us perform Fourier transformation with 
respect to the spatial coordinate r. The concentration after both transformations is denoted 
by c(k, z). Making the ansatz c(r, t) = g(r) + <p(r, t) or 

c(k, z) = ^ + ^(k, z) (7) 

z 

where <p(\&, z) = C{ip(k, t)} = J °° </?(k, £)} exp(— zt)dt should not display a singular behavior 
for z —>■ 0. Inserting Eq. (0) into Eq.(0]) the singular part yields the stationary solution c s (k) 
which obeys the equation 

c s (k)[(X-^k 2 )c (k)]-ub(k) 
° s[ ' Dk 2 + {\- /iP)c s (k) 

with 6(k) = y^c,(q)c s (k-q), c„(k) = c(k,t = 0) (8) 

This relation has always the solution c s (k) = which corresponds to a complete decay, 
ending up with zero concentration. In addition to this trivial solution Eq. (jHJ) exhibits 
a non-zero stationary solution, although the instantaneous decay process, manifested the 
decay rate u 7^ 0, is still present. Obviously the memory effects prevent a complete reaction. 
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The non-trivial solution of Eq. (JHJ) depends on the initial value co(k) and on the parameters 
of the model. A linear stability analysis in terms of <£>(k, t), Eq. (It), leads to 

V?(k,t) = -A(k)p(k,*)-2u [ c s {k - q)<p{q,t) 

Jq 

with A(k) = Dk 2 + (A - fik 2 )c s (k) (9) 

To discuss the generic behavior let us firstly study the case u = which is reasonable, 
since the inhomogeneous stationary solution is originated mainly by the memory terms 
proportional to the coupling parameters \i and A, respectively. In that case the nontrivial 
stationary solution reads in according to Eq. (jHJ) 

c s (k; u = 0) = c (k) - (10) 

If A = it results a homogeneous solution as expected whenever the initial concentration is 
likewise homogeneous. In general the stationary solution is only accessible when c s (k) > 
and reasonable if c s (k; u = 0) < c (k). Because the asymptotic behavior depends essentially 
on the signs of the feedback parameters fi and A we consider firstly the case of fi > and 
A > 0. For positive parameters \i and A > 0, the inhomogeneous stationary solution is stable 
and reasonable for k < k c where the critical wave vector is determined by 

2 _ Acq (fe c ) 

fcc D + ^c (k c ) [LL) 

A simple realization is given for c(r, t = 0) = c <5(r) leading to c (k) = c . In that case 
the critical wave vector is simply defined by the last equation where c(k c , 0) is replaced by 
Cq. Remark that the critical value k c fulfills the relation k c < \J\j \i In this manner the 
spurious singularity in Eq. ()10|) is avoided. Together with the stability criteria A(k) > we 
get the phase diagram depicted in Fig^ The curve is similar to the behavior of a second 
order phase transition. Introducing the variable w = k/k c we get 



, c (£> + c /i) 

^;n = 0) = (i- w ) D + cMi _ w ^ (12) 



One can study the spatial dependence of non-trivial behavior of c s (r) by making the inverse 
Fourier-transformation. For large r =| r |— > oo we get 

c s (r) = c(r,0) 



d+2 
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where k is a factor depending on the parameters of the model. The physical behavior of 
the system can be discussed in terms of particles and holes denoted as A and B particles, 
respectively. On a large length scale the final state of the system consists of homogeneous 
distributed A-particles, which becomes inert and which are separated from the holes. On a 
short scale the A-particles disappear further until they are completely annihilated and only 
holes leave over. The systems exhibit a new kind of phase separation mechanism. 
In the same manner we can discuss the residual cases with varying signs of the feedback 
parameters /i and A. If both parameters are negative, fj, < and A < 0, we find the inverse 
situation compared to that one discussed above. A non-zero stable and simultaneously 
reasonable stationary state appears for k > k c , whereas for k < k c all the particles undergo 
the complete reaction, see FigO The system exhibits a microphase-separation. It decays 
in small subunits of extension I < k~ x . If /j > and A < there exists only the trivial 
stationary solution c s (k) = 0. In the remaining case, /x < and A > 0, we find a stable 
non-zero stationary solution for all wave vectors provided the initial concentration fulfills 
the condition co(k) > D | /i In the opposite case one observes a similar behavior as it is 
depicted in FigHJ 

If the reaction parameter u is non-zero, Eq. (jHJ) can not be solved exactly due to the non- 
linear convolution term. However, the main conclusions remain valid, in particular for a 
small coupling parameter u. In a linear approximation with respect to u the stationary 
solution Eq. (fT0|) is modified in according to c s (k; u = 0) — > c s (k; u) = c s (k; u — 0) — uR(k) 
with 

= / d d q/(2n) d ) [ Ca (q; u = 0)c s (k - q; u = 0)] 
1 ' (A — fxk 2 )c s (k; u = 0) 

Expanding the last expression with respect to k 2 by using Eq. (fTU|) we find, that the critical 
wave vector k c , compare Eq. (fTTj). is slightly reduced and the stability of the solution is 
maintained. The phase diagram, presented in Fig^ is preserved. 

IV. CONCLUSIONS 

In this paper we have extended the conventional modeling of chemical reactions by includ- 
ing non-Markovian memory terms within the evolution equation. The additional terms give 
rise to a competitive behavior for both the local reaction and the diffusive transport term. 
The influence of the feedback-couplings is essential, in particular in the long time limit. Due 



S 



to the memory effects the reactants may be localized and therefore a further reaction is pre- 
vented. The resulting asymptotic inert state is stable on a large length scale. However, there 
exists a critical wave vector, above it, the reaction can be continued as long as all the react- 
ing particles are annihilated completely . The reason for such a new behavior is by means 
of an explicit coupling of the rate of the concentration at the observation time t to that one 
at a previous time. This time accumulation is further accompanied by an additional spatial 
accumulation, the effect of which is comparable to the effect a long-range interaction forces 
and consequently the results are basically independent on spatial dimensions in according to 
scaling arguments. These many-body effects are shown to change the asymptotic behavior 
drastically. Due to the feedback-coupling of a particle to its environment, a subsequent par- 
ticle, undergoing a diffusive motion, gains information from a modified environment. Hence 
the particle can be confined within a certain region preventing a further reaction. In this 
manner a self-organized memory leads to a non-zero stationary concentration c s (r), where 
the analytical form of c s is controlled by the memory strength. This situation is realized on 
a large length scale. Below a critical wave vector k c the inhomogeneous stationary solution 
is stable. The phase diagram is similar to an equilibrium phase transition of second kind. 
The model could be relevant for very complex chemical reactions where not all entities are 
present simultaneously. 
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FIG. 1: Stationary concentration c s (k) versus w = k/k c for positve parameters fi, A > 
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FIG. 2: Stationary concentration c s (k) versus w = k/k c for negative parameters < and A < 
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